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Abstract 

This paper presents explicit algebraic transformations of some Gauss hypergeo- 
metric functions. Specifically, the transformations considered apply to hypergeomet- 
ric solutions of hypergeometric differential equations with the local exponent differ- 
ences 1/k, l/£, 1/m such that k, £, m are positive integers and 1/k + 1/i + 1/m < 1. 
All algebraic transformations of these Gauss hypergeometric functions are consid- 
ered. We show that apart from classical transformations of degree 2, 3, 4, 6 there are 
several other transformations of degree 6, 8, 9, 10, 12, 18, 24. Besides, we present 
an algorithm to compute relevant Belyi functions explicitly. 
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1 Introduction 

An algebraic transformation of Gauss hypergeometric functions is an identity of the form 
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where (p(x) is a rational function of x, and 9(x) is a product of some powers of rational 
functions. Here are two examples of quadratic transformations (see [Erd53, AAR99]): 
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These identities hold in some neighborhood of x = in the complex plane, and can be 
continued analytically. For example, formula (2) holds for Re(a;) < 1/2. 
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Recall that the Gauss hypergeometric function 2 Fi 
pergeometric differential equation 




is a solution of the hy- 
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This is a Fuchsian equation on the complex projective line P 1 with 3 regular singular 
points z = 0, 1 and oo. The local exponent differences at these points are (up to a sign) 
1 — C, C — A — B and A — B respectively. 

Algebraic transformations of Gauss hypergeometric functions usually come from those 
transformations of hypergeometric equation (4) , which have the form 



and such that the transformed equation for Y{x) is a hypergeometric equation in the new 
indeterminate x. Here <p(x) and 8(x) have the same meaning as in formula (1). Geometri- 
cally, this is a pull-back transformation of equation (4) with respect to the finite covering 
if : P 1 — > P 1 determined by the rational function <p(x). In [Kit03] these transformations 
are called RS-transformations. Recall that a rational function on a Riemann surface is 
a Belyi function [ShaOO, Krc03] if it has at most 3 critical values, or equivalently, if the 
corresponding covering of P 1 branches only above a set of 3 points. The function ip(x) in 
hypergeometric identities like (1) is usually a Belyi function. 

Algebraic transformations of Gauss hypergeometric functions and pull-back transfor- 
mations of hypergeometric equations are related as follows. 

Lemma 1.1 1. Suppose that pull-back transformation (5) of equation (4) is a hyper- 
geometric equation as well, and that the transformed equation has non-trivial mon- 
odromy. Then, possibly after fractional-linear transformations on the projective 
lines, there is an identity of the form (1) between hypergeometric solutions of the 
two hypergeometric equations. 

2. Suppose that identity (1) holds in some region of the complex plane. Let Y(x) 
denote the left-hand side of the identity. If Y 1 (x) /Y{x) is not a rational function 
of x, then the transformation (5) converts the hypergeometric equation (4) into a 
hypergeometric equation for Y (x) . 

Proof. This is Lemma 2.1 in [Vid04]. □ 

In this paper, we consider Gauss hypergeometric functions which satisfy hypergeomet- 
ric equations with local exponent differences 1/fc, l/l, 1/m such that k,£,m are positive 
integers and 1/k + l/i+l/m < 1. We call these functions hyperbolic hypergeometric func- 
tions, because they have interesting analytic properties related to the hyperbolic geometry 
of the complex plane [Yos97, Beu02] . The main purpose of this paper is to describe alge- 
braic transformations of these functions into other hypergeometric functions. Existence 
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of their non-classical transformations of degree 10, 12 and 24 is shown in [Hod20, Bcu02]. 
A transformation of degree 8 is presented in [Kit03, Section 5]. We give a complete list 
of possible algebraic transformations of hyperbolic hypergeometric functions. Algebraic 
transformations of all Gauss hypergeometric functions are classified in [Vid04] . 

For hyperbolic hypergeometric functions, algebraic transformations always induce 
pull-back transformations of their hypergeometric equations, and vice versa. Indeed, Ko- 
vacic algorithm [Kov86, vdP98] in differential Galois theory implies that the monodromy 
group of those hypergeometric equations is not trivial, and that they have no solutions 
y(z) with algebraic logarithmic derivative y'(z)/y(z). Therefore Lemma 1.1 allows no 
exceptions. 

This paper classifies algebraic transformations of hyperbolic hypergeometric functions 
by finding all pull-back transformations of their hypergeometric equations to other hy- 
pergeometric equations. The main problem is to compute suitable coverings ip : P 1 — > P 1 . 
Possible branching patterns for them are derived in Section 2. We give a general algo- 
rithm for computing coverings with prescribed branching pattern in Section 3. Algebraic 
transformations of hyperbolic hypergeometric functions are listed in Section 4. 

2 Possible branching patterns 

A general pull-back transformation (5) of a hypergeometric equation is a Fuchsian equa- 
tion. We are looking for situations when the transformed equation is hypergeometric 
as well. In this Section we rather look for transformed equations with at most 3 singu- 
lar points. Since any such Fuchsian equation can be transformed to a hypergeometric 
equation by fractional- linear transformations, it is appropriate to ignore exact location of 
singular points for a while. We loosely follow the 5-step classification scheme in [Vid04, 
Section 3], with N = 3, etc. 

The requirement that the transformed equation must have at most 3 singular points 
is restrictive. The covering ip : P 1 — > P 1 essentially determines singularities and local 
exponent differences of the transformed equation. Here are basic general facts which we 
use (or refer to). 

Lemma 2.1 Let ip : P* — > Pj denote a finite covering of a projective line with the 
rational parameter z by a projective line P* with the rational parameter x. Let Hi denote 
a hypergeometric equation on ¥],, and let Hi denote the pull-back transformation of Hi 
under (5). Let d denote the degree of ip, and let S E P* ; Q € P* be points such that 

v(S) = Q- 

1. If the point Q is non-singular for Hi, then the point S is non- singular for H 2 only 
if the covering ip does not branch at S. 

2. If the point Q is a singular point for Hi , then the point S is non-singular for Hi 
only if the local exponent difference at Q is equal to l/n, where n is the branching 
order of ip at S. 
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3. Let A denote a set of 3 points on P*. If all branching points of ip lie above A, then 
there are exactly d + 2 distinct points on above A. Otherwise there are more 
than d + 2 distinct points above A. 

4- Suppose that the equations Hi and H 2 are hypergeometric. Let e\,e 2 ,e3 denote 
the local exponent differences for H\, and let e'i,e' 2 ,e' 3 denote the local exponent 
differences for H 2 . Suppose that the local exponent differences are real positive 
numbers, and that e\ + e 2 + e 3 ^ 1. Then 



1 



ei - e 2 



e3 



(6) 



Proof. The first two statements are weaker formulations of parts 2, 3 of [Vid04, Lemma 
2.4]. The third statement is part 1 of [Vid04, Lemma 2.5], and the last statement is a 
weaker formulation of part 2 of [Vid04, Lemma 2.5]; they are consequences of Hurwitz' 
formula [Har77, Corollary IV.2.4]. □ 

Here are restrictions on coverings p> : P 1 — > P 1 and local exponent differences for alge- 
braic transformations of hyperbolic hypergeometric functions. They are stronger versions 
of the constraints used in Step 2 of the classification scheme in [Vid04, Section 3] . 



Lemma 2.2 Let k,£,m denote positive integers such that 

and k < i < m. 
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Let H denote hypergeometric equation (4) such that the local exponent differences are 
equal to l/k,l/£,l/m. Suppose that pull-back transformation (5) transforms H to a 
hypergeometric equation. Let d denote the degree of the covering p : P 1 — > P 1 . 

1. The points x = 0, 1, oo are actual singularities of the transformed equation, and they 
lie above the subset {0, 1, oo} of the z-projective line. The covering ip branches only 
above this subset, so <p(x) is a Belyi function. 



2. The following equality holds: 
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3. The following inequality holds: 
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I If m> d then 1/d+l/k +1/£>1. 
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5. If m < d then 
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Proof. Let A denote the subset {0, 1, oo} of the ^-projective line. By part 2 of Lemma 
2.1, there are at most [d/k\, [d/£\, [d/m\ non-singular points above the z-points with 
the local exponent differences 1/k, l/£, 1/m respectively. By part 3 of Lemma 2.1, the 
number of singular points above A is at least 
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This number is greater than 2 + rf(l — 1/k — l/£ — 1/m), so it is at least 3. On the other 
hand, the transformed equation has at most 3 singular points. Therefore the transformed 
equation has exactly 3 singular points, expression (11) is equal to 3, and ip does not 
branch outside A by part 1 of Lemma 2.1. Parts 1 and 2 of this Lemma follow. 
We rewrite formula (8) as follows: 



d 1 - 



1 



+ T = 1, 



(12) 



where T is the sum of positive local exponent differences at the singular points of the 
transformed equation. (This is equivalent to (6), with T — e[ + e' 2 + e' 3 .) We have 
T > 3/m, which implies inequality (9). 

If m > d then we use formula (8) to derive 



(13) 



which gives part 4 of this Lemma. If m > d, we derive the first inequality in (10) after 
replacing d by m in (9). We have l/k+l/£ < 1 from (7). The quadratic expression in m 
achieves non-positive values only ifl/fc + l/^>2/3 (consider the discriminant). □ 
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These restrictions essentially give a finite list of possibilities for the integer tuple 
(k,£,m,d). Indeed, inequality (9) bounds d once k,£,m are fixed, and part 5 of Lemma 
2.2 gives finitely many possibilities for the triple (k,£,m). Only when m > d we formally 
have infinitely many possibilities; but then we expect to arrive at specializations of alge- 
braic transformations with unrestricted parameters. The inequalities give the following 
possibilities: 

(2,£,m,2), (2,3,m,3...6), (2,4,m,4), (3,3,m,3), (2, 3, 7, 7...24), 
(2,3,8,8...15), (2,3,9,9...12), (2,3,10,10), (2, 4, 5, 5...S), (2,4,6,6). 

Here d is sometimes represented by an integer interval of possible values, and the uneval- 
uated parameters £, m can be large enough integers. Formula (8) rejects some of these 
possibilities. 
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Local exponent differences 
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Covering 
composition 


Coxeter 
decomposition 


C1 Ik 1 IP 1 Im) 


au\j v vi- 


(1/2, 1/*, 1/m) 


(1/4 1/4 2/m) 


Z 


indecomposable 


yes 


(1/2, 1/3, 1/m) 


(1/2, 1/m, 2/m) 


3 


indecomposable 


yes 


(1/2, 1/3, 1/m) 


(1/3, 1/m, 3/m) 


4 


indecomposable 


yes 


(1/2, 1/3, 1/m) 


(1/3, 2/m, 2/m) 


4 


no covering 


(1/2, 1/3, 1/m) 


(1/m, 1/m, 4/m) 


6 


2x3 


yes 


(1/2, 1/3, 1/m) 


(2/m, 2/m, 2/m) 


6 


2 x 3 or 3 x 2 


yes 


(1/2, 1/3, 1/m) 


(1/m, 2/m, 3/m) 


6 


no covering 


(1/2, 1/4, 1/m) 


(1/m, 1/m, 2/m) 


4 


2x2 


yes 


(1/3, 1/3, 1/m) 


(1/m, 1/m, 1/m) 


3 


indecomposable 


no 



Table 1: Classical transformations of hyperbolic hypergeometric functions 



The next step is to produce a list of possible branching patterns. Because of parts 1 
and 2 of Lemma 2.2, we have to take the maximal possible number [d/k\, [d/P\ or \d/m\ 
of non-singular points above the 3 singular ^-points. The remaining residual branches 
above z = 0, 1, oo should coalesce into precisely 3 distinct points. In particular, we have 
to ignore the cases when there remains less than 3 residual branches. The final list of 
possible branching patterns is presented in the first three columns of Table 1 (for the cases 
with m > d) and Table 2 (for the cases with m < d). In Table 1 we ignore fractional- 
linear transformations, and we drop the condition P < m for degree 2 transformations. 
Branching patterns are uniquely determined by the starting local exponent differences 
(l/k,l/P,l/m), transformed local exponent differences, and the stated principle to have 
maximal number of non-singular points above z = 0, 1, oo. For example, branching pat- 
tern for the degree 4 transformation of Table 1 between hypergeometric equations with 
the local exponent differences (1/2, 1/3, 1/m) and (1/3, 1/m, 3/m) can be schemati- 
cally denoted by 2 + 2 = 3 + 1 = 3 + 1. This means that all points above the z-point 
with the local exponent difference 1/2 have branching order 2, and that there must be 
a branching point with order 3 and a non-branching point above each of the other two 
points. For one more example, the degree 9 covering of Table 2 has the branching pattern 
2 + 2 + 2 + 2 + 1 = 3 + 3 + 3 = 7+1 + 1. The same notation for branching pattern is 
used in [Vid04]. 

Now we have to determine all coverings which have those branching patterns. Given 
a branching pattern, there is often exactly one covering with that branching pattern up 
to fractional-linear transformations. But not for any branching pattern a covering exists, 
and there can be several different coverings with the same branching pattern. Section 3 
is devoted to computing coverings with a given branching pattern. First we outline there 
a straightforward method with undetermined coefficients, which is feasible if d < 6. Then 
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Local exponent differences 
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1 A 
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It It i It o It\ 

(1/7, 1/7, 3/7) 




no coverin; 
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(1/2, 1/3, 1/7) 


(1/7, 2/7, 2/7) 




no coverin; 
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(1/2, 1/3, 1/7) 


(1/3, 1/7, 1/7) 


16 


no coverin; 
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(1/2, 1/3, 1/7) 


(1/7, 1/7, 2/7) 


18 


2x9 




no 


(1/2, 1/3, 1/7) 


(1/7, 1/7, 1/7) 


24 


3x8 




yes 


(1/2, 1/3, 1/8) 


(1/3, 1/8, 1/8) 


10 


indecomposable 




no 


(1/2, 1/3, 1/8) 


(1/4, 1/8, 1/8) 


12 


2x2x3 




yes 


(1/2, 1/3, 1/9) 


(1/9, 1/9, 1/9) 


12 


3x4 




no 


(1/2, 1/4, 1/5) 


(1/4, 1/4, 1/5) 


6 


indecomposable 




no 


(1/2, 1/4, 1/5) 


(1/5, 1/5, 1/5) 


8 


no coverin; 


r 



Table 2: Non-classical transformations of hyperbolic hypergeometric functions 



we introduce a more appropriate algorithm, which was actually used (within computer 
algebra system Maple) to compute coverings for Table 2. Final information about existing 
coverings is given in the fourth columns of Table 1 and Table 2. It turns out that for 
any candidate branching pattern there is at most one covering up to fractional-linear 
transformations. If m > d (see Table 1), we get the coverings of the classical algebraic 
transformations due to Gauss, Euler, Kummer, Goursat. If m < d (see Table 2), we 
get new coverings of degree 6, 8, 9, 10, 12, 18, 24. Existence of some of these coverings 
is shown in [Hod20, Bcu02, Kit03]. For both Tables, it was straightforward to figure 
out possible compositions of small degree coverings and identify them with the unique 
coverings for suitable branching patterns. Numbers in the multiplicative notation for 
decomposable coverings mean degrees of constituent coverings, as in [Vid04] . 

The last step is to determine algebraic transformations of hypergeometric functions 
with the rational argument determined by a computed covering. The factor 9(x) in (5) 
should shift local exponents at potentially non-singular points to the values and 1, and 
it should shift one local exponent at both x = and x = 1 to the value 0. A suitable pull- 
back transformation induces a hypergeometric identity like (1) for each singular x-point 
S which lies above a singular z-point. To achieve this, one has to move the points S and 
(p(S) to the locations x = and z = respectively (by fractional- linear transformations), 
and identify the two solutions with the local exponent and the value 1 at x = and 
z = respectively. It is convenient to use Ricmann's P-notation for these purposes; see 
[AAR99, Section 3.9] or [Vid04, Section 2]. Each positioning of x — above z = gives a 
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few hypcrgcomctric identities like (1). First of all, we have Euler's and Pfaff's fractional- 
linear transformations [AAR99, Theorem 2.2.5], which permute other two singular points 
and their local exponents. Additionally, simultaneous permutation of the local exponents 
at x = and z = gives the following hypergeometric identity. 

Lemma 2.3 Suppose that a pull-back transformation induces identity (1) in an open 



neighborhood of x = 0. Then if(x) 1 c ~ Kx 



,i-C 



as x — » for some constant K , and the 



following identity holds (if both hypergeometric functions are well-defined): 
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Kx 



1-C 



2F1 



1 + A-C,1 + B-C 

2-C 



p(x] 



Proof. This is Lemma 2.3 in [Vid04]. 



□ 



As it turns out, algebraic transformations for Table 1 (i.e., the case m > d) are special 
cases of the classical transformations due to Gauss, Euler, Kummer, Goursat. We give a 
few instances of these transformations in Section 4. Algebraic transformations for Table 
2 (i.e., the case m < d) arc modern, though some of them are predicted in [Hod20]. We 
present these transformations (up to Euler's and Pfaff's fractional-linear transformations, 
and Lemma 2.3) in Section 4 as well. 

The rest of this Section is devoted to explaining the last columns of Tables 1 and 
2. Recall [Yos97, Beu02] that a Schwarz map for a hypergeometric equation H\ is an 
analytic map from the upper half-plane H = {z G C | Im z > 0} given by a quotient of 
two solutions of Hi. If the local exponent differences eo,ei,eoo of Hi are real numbers 
in the interval [0, 1), then the image of a Schwarz map is a curvilinear triangle on the 
Ricmann sphere. Such a triangle is called Schwarz triangle. The vertices are images of 
the 3 singular points, and the angles there are equal to ne , irei, -ne^ correspondingly; 
the sides are circular arcs. Analytic continuation of a Schwarz map follows the Schwarz 
reflection principle: the image of the other half-plane under analytic continuation across 
(0, 1), (1, oo ) or (—00, 0) is a fractional-linear reflection of the Schwarz triangle across the 
corresponding side of itself. 

In our case, the local exponent differences are 1/fc, l/£,l/m, and 1/k + l/l+l/m < 1. 
The sides of a Schwarz triangle are geodesic curves with respect to a hyperbolic metric on 
the Riemann sphere, defined on some Poincare disk. Repeated analytic continuation gives 
a tessellation of the Poincare disk into curvilinear triangles with the angles n/k, n/l^/m. 

Consider a pull-back transformation of the hypergeometric equation Hi to a hyperge- 
ometric equation H 2 , of degree d. Suppose that its covering ip : P 1 — ► P 1 is defined over 
K, and that it branches only above the singular points of Hi . If s : HI — > C is a Schwarz 
map for H 2 , then a branch of s o ip^ 1 is a Schwarz map for H±. The Schwarz triangles of 
the d branches of s o ip' 1 tessellate the Schwarz triangle of s, like in Figure 1. In this case 
the degree expression (6) can be interpreted as the quotient of areas of Schwarz triangles 
for the two hypergeometric equations, in the hyperbolic or spherical metric. Transforma- 
tions of hypergeometric equations which admit these tessellations are implicitly classified 
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Figure 1 : Coxeter decompositions of hyperbolic triangles 



in [Hod20, Beu02]. Tessellations of hyperbolic triangles and quadrangles into hyperbolic 
triangles are classified in [Fel98, BHMM], where they are called Coxeter decompositions 
and divisible tilings respectively. The classification in [BHMM] is incomplete; for instance, 
it misses triangulation (b) in Figure 1. We adopt the terminology of [Fel98]. 

The last columns of Tables 1 and 2 tell us which transformations of hypergeometric 
equations with hyperbolic solutions admit Coxeter decompositions of Schwarz triangles. 
In particular, the three such transformations in Table 2 are anticipated in [Hod20, Bcu02]. 
Their Coxeter decompositions are depicted in Figure 1. All classical transformations 
except one cubic transformation admit these tessellations. 

3 Computation of Belyi functions 

Here we consider the problem of computing finite coverings ip : P 1 — > P 1 of given degree 
d and with a given branching pattern. We assume that all branching points lie above 
a set A of 3 points, so the desired functions f(x) are Belyi functions. By part 3 of 
Lemma 2.1, there must be exactly d+2 distinct points above A. We expect finitely many 
(or no) solutions to this problem. Algorithm 3.1 of this Section was used to compute 
transformations implied by Tables 1 and 2. As mentioned, there is at most one solution 
for branching patterns there. 
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ya 2 j, T «1A t «o; / 1 t \ 

^ = ^27— n2 . ( 14 ) 



First we outline a naive method with undetermined coefficients, which is feasible if 
d < 6. To fix ideas, consider the branching pattern 2 + 2 + 2 = 3 + 3 = 2 + 2 + 2 for 
the sixth entry in Table 1. Up to fractional- linear transformations of the 2-line, we may 
assume that points with these branching orders lie above z = 1,0 and oo respectively. 
We choose the points above z = oo to be x = 0, x = 1, x = oo. Then the Belyi function 
should have the form 

(u 2 x 2 + u\X + u ) 3 

x 2 (x — iy 

where u 2 , Mi, uo are undetermined, and the roots of u 2 x 2 + Uix + uq are the points above 
z = 0. The branching pattern above z = 1 implies that the numerator of ip{x) — 1 must 
be a square of a cubic polynomial P{x). This condition gives 7 polynomial equations 
in the 4 coefficients of P(x) and 112,1*1, wo- These equations can be feasibly solved with 
assistance of a computer algebra package if d is not large. 

To compute Belyi functions more efficiently, we propose to pull-back the differential 
dz/z with respect to <p, still with undetermined coefficients as in (14). The poles of the 
pull-backed differential are simple, and they are located at the points above z = and 
z = 00. Its zeroes are the branching points which do not lie above z = or z = 00; 
the multiplicities of the zeroes are 1 less than the corresponding branching orders. In 
our example, all those branching points must lie above z = 1, so they are the roots of 
P(x). Moreover, they must be simple roots of P(x), so we get this polynomial just by 
computing the pull-back of dz/z. Explicitly, the pull-back of dz/z is equal to 

ifi'jx) _ 2u 2 x 3 - (4m 2 +mi) x 2 - (mi+4m ) x + 2u 

; — CLX — - . . 7. T ClX . ( 1 j 

<p{X) X (X — 1) (ll2X z + U\X + Uo) 

Let P(x) denote the polynomial in the numerator on the right-hand side; it must be 
proportional to P{x). Therefore <p(x) — 1 is proportional to <p(x) :— P(x) 2 / x 2 (x — l) 2 . 
Further, consider the pull-back of the differential d(z — l)/(z — 1): 

<p'(x) _ , u 2 x 4 — 2u 2 x 3 + (2u 2 +ui + 2u )x 2 — 2uqx + u 
ip(x) x (x — 1) (2u 2 x 3 ~ (4:U 2 + ui)x 2 — (ui+Auo) x + 2uo) 

By the same reasoning, the zeroes of this differential are the branching points of ip which 
do not lie above z = 1 or z = 00, with the multiplicities diminished by 1. In our case, 
those branching points must lie above z = 0. Hence the polynomial in the numerator of 
(16) is proportional to (u 2 x 2 + u\x + uq) 2 . This gives easy polynomial equations in u 2l 
u\, uq. Since we want 1*2^0 7^ 0, 

u\ _ 2u x u 2 _ u\ + 2u Q u 2 _ 2w ui _ u 2 , 
u 2 2u 2 2u 2 +u\ + 2uq 2uq uq 

We solve that u 2 = —u\ = Mo- Therefore ip(x) is proportional to (x 2 — x + l) 3 /a; 2 (.x — l) 2 . 
The scalar multiple can be found from the condition that z = 1 is the third branching 
locus of ip(x). We derive: 

^ > 27 x 2 (x~\) 2 ' V( ' 27x 2 (x-l) 2 
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We have solved the problem by hand! The solution is unique up to fractional-linear 
transformations. Note that there are two different ways to compose coverings of degree 2 
and 3 and get a covering with the considered branching pattern; see Table 1 and [Vid04, 
Section 4] . Up to fractional- linear transformations, those two compositions must give the 
same covering computed here. 

Now we present general Algorithm 3.1 for finding Bclyi functions with a given branch- 
ing pattern. To formulate it more conveniently, we restrict ourselves to branching patterns 
that are relevant for the purposes of this paper. Note that for all transformations for Table 
2 (and almost all transformations for Table 1) there is a z-point with the local exponent 
difference 1/2. For coverings of these transformations assumption (b) of Algorithm 3.1 
holds. If this assumption is dropped, then Step 1 should try to assign the fiber with 
smallest branching orders to z = 1, the function ip(x) in Step 3 has a more complicated 
form, and more undetermined coefficients are needed. 

Algorithm 3.1 Input: a branching pattern (that is, 3 collections of branching orders) 
and degree d. We assume: 

(a) the branching orders in the same fiber sum up to d, and there are d + 2 branching 
orders in total; 

(b ) one of the 3 collections prescribes only branching orders 2 and at most one unramificd 
(i.e., simple, not branching) point. 

Output: All Belyi functions (up to fractional-linear transformations) whose coverings 
branch only above the set A = {0, 1, oo} with the given branching orders. 
Step 1. Prescribe the branching orders mentioned in assumption (b) to the fiber of 
the point z = 1, and prescribe other two collections to the fibers of z — and z = oo. 
Choose the points £ = 0, x = 1, x = oo above A in a convenient way: if an unramificd 
point is prescribed above z = 1, choose it to be x = oo; see also remarks immediately 
below. Consider the other points above z = and z = oo as unknown. Accordingly, 
write ip(x) = K P(x)/Q(x), where K is an undetermined constant, and P(x), Q(x) are 
monic polynomials in the square-free factorized form (following the branching pattern) 
with some undetermined coefficients. 

Step 2. Compute the pull-back ip'(x) dx/ip(x) of dz/z. Let R(x) be the numerator of 

if' (x) / tp(x) . The roots of R(x) are the branching points above z = 1. 

Step 3. Let tp(x) = R(x) 2 /Q(x) and compute the rational function &(x) = tp'(x)/ip(x). 

The numerator of $(x) has the same roots as the polynomial P(x), but their multiplicity 

is 1 less than in P(x). This gives a set of algebraic equations in the undetermined 

coefficients. 

Step 4. Solve the algebraic equations obtained in Step 3 by Grobner basis methods, 
and find possible pairs of polynomials P(x),Q(x) with the right factorization pattern. 
For each non-degenerate solution, the constant K in the target (p(x) = K P(x)/Q(x) is 
such that the function 1 — K P(x)/Q(x) is proportional to <f(x). If necessary, compose 
the output functions <p(x) with suitable fractional-linear transformations to move some 
appoints (or even ^-points) to final desired locations. □ 
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When applying this algorithm to the entries of Table 2, it is convenient to choose the 
points x = 0, x = 1, x = oo in Step 1 to be the singular points of the transformed 
hypergeometric equation. In general, a good strategy for Step 1 is to choose points which 
have different branching orders than the most points in the same fiber. On the other 
hand, the algorithm can be more effective if we choose points with maximal branching 
orders as x = 0, x = 1, x = oo. With this modification, the function tp(x) in Step 3 
may acquire an extra linear factor. It may be convenient not to make a choice for x = 1 . 
Then we would have an extra variable and the algebraic equations would be weighted- 
homogeneous (respecting to the transformations x i— > ax). The extra degree of freedom 
can be used to avoid complicated algebraic numbers. 

We note that Step 3 produces enough algebraic equations between the undetermined 
coefficients, because the restrictions on the polynomials P(x), Q(x), R(x) and the denom- 
inator of $(ar) determine the desired branching pattern for <p(x). Hence the algorithm is 
correct. As example (17) shows, the set of equations is likely to be overdetermined, which 
only helps in Grobner basis computations. 

Compared with the naive method described at the beginning of this Section, algebraic 
equations of Algorithm 3.1 have fewer undetermined coefficients, lower degree, and fewer 
degenerate (or parasytic [Krc03]) solutions. The computations are still tedious, but all 
coverings of Table 2 were computed using the computer algebra package Maple in a matter 
of hours. There is an article [ShaOO] where quadratic differentials are used to characterize 
some Belyi maps. But [Kre03, Kit03] exploit the naive method. 



4 Hypergeometric identities 

Here we present our main results. We give all algebraic transformations for Table 2, up to 
Eulcr's and Pfaff's fractional-linear transformations and Lemma 2.3. But first we exhibit 
a few classical transformations. 

Relevant instances of quadratic transformations can be obtained by setting a = 1/2 — 
l/£ — l/m,b= 1/2 — 1/£ +l/m in formulas (2)-(3). Examples of classical transformations 
of degree 3 or 4 are: 



2F1 
2 Fi( 



I ia+5 
\ 6 

4a Aa+l 
3 ' 3 
4a+5 



2F1 



- £+1 

3 

2c+2 



X = 



(1 - 4xy a 2 Fi y 
(l + 8xy a 2 F^ 
(l+uj 2 X y c 2 Fi| 



a a+1 
3 ' 3 


27 a; \ 


4a+5 
6 


{4x - i) 3 y ' 


a a+1 
3 ' 3 


64 a: (1 - xf 


4a+5 
6 


(1 + 8x) 3 



c e+1 
3' 3 
2c+2 
3 



3(2w + l) x(x - 1) 
{x + w) 3 



(18) 
(19) 
(20) 



Here w is a primitive cubic root of unity (so uo 2 + u> + 1 = 0), a = 1/4 ± 3/2m and 
c = 1/2 ± 3/2m. These identities correspond to the indecomposable pull-back coverings 
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if : P 1 — > P 1 of Table 1 . For more complete lists of classical algebraic transformations we 
refer to [Gou81, Vid04]. 

Now we present non-classical transformations of hyperbolic hypergeometric functions. 
Our list of transformations is basically complete, it was computed following the plan 
of Section 2. If a covering for Table 2 is indecomposable, all corresponding two-term 
hypergeometric identities can be obtained from the exhibited below, by using Euler's and 
Pfaff 's fractional-linear transformations and Lemma 2.3. If the covering is decomposable, 
we indicate a composition of hypergeometric identities of smaller degree. 

A covering for degree 8 pull-back transformations between hypergeometric equations 
with the local exponent differences (1/2, 1/3, 1/7) and (1/3, 1/3, 1/7) is given by 



ipx{x) = 



x (x - 1) (27a; 2 - (723 + 1392w)a; - 496 + 696w) 3 
64((6w + 3)a;-8-3w) 7 



(21) 



Here ui satisfies uj 2 + lo + 1 = as in formula (20). Note that the conjugation u = — 1 — u> 
acts in the same way as a composition with fractional-linear transformation interchanging 
the points x — and x = 1. This confirms uniqueness of the covering. The covering is 
computed in [Kit03] as well. Here are hypergeometric identities: 



2F1 
2F1 



2/21, 5/21 

2/3 
/2/21, 3/7 
I 6/7 



(1 

(1- 

X 2 F! 



33+39^-1/12^ 



49 



1/84, 13/84 
2/3 



241+464^ _ 8(62-87^) 2 
9 27 



) 



■1/28 



/ 1/84, 29/84 
I 6/7 



1 



ipi(l/x) 



(22) 



(23) 



A covering for degree 9 pull-back transformations between hypergeometric equations 
with the local exponent differences (1/2, 1/3, 1/7) and (1/2, 1/7, 1/7) is given by 



ip 2 {x) 



27x(a;-l)(49a;-31-130 7 



49 (7203a; 3 + (9947£ - 5831)a; 2 - (9947£ + 2009)a; + 275 - 87£) 3 ' 



(24) 



Here £ satisfies £ 2 + £ + 2 = 0. Hypergeometric identities are: 



2F1 



/3/28, 17/28 
I 6/7 



2F1 



( 3/28, 1/4 
I 1/2 



7(10-29g) 



X — 



343(50-29g) 



X = 



X2F1 
X2F1 



/ 1/84, 29/84 
I 6/7 



512 



x' + 



1029(362+87g) 3 



-1/28 



16384 



X 



17-294 _ 41+203$ 2 , 275-87g 3 
21 ■*' 147 ~ 7203 , 



-1/28 



1/84, 29/84 
1/2 



l-V2(l/aO 



(25) 



(26) 
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A covering for degree 10 pull-back transformations between hypergeometric equations 
with the local exponent differences (1/2, 1/3, 1/7) and (1/3, 1/7, 2/7) is given by 



x 2 (x- l)(49a;-81) 7 



4 (16807a; 3 - 9261a; 2 - 13851a; + 6561) 3 



Hypergeometric identities are: 
2F1 



2F1 



2*1 



5/42, 19/42 
5/7 
/5/42, 19/42 
I 6/7 



1 19 343 2 

x g a. 243 ^. 



16807 ^,3\~ 1 / 28 p 
~&XT X ) 2 *1 



a; 



(l-Iflz+^a 



16807 
256 



X 3 ) 



-1/28 



Fi 



1/84, 29/84 

6/7 
/ 1/84, 29/84 
I 6/7 



(27) 



V3(a0 , (28) 



/5/42, 17/42 
I 2/3 



{1 - x )-^(l-§xy 1/12 2 F^ 



( 1/84, 13/84 
2/3 



^3(1/2;) 



ip 3 (l-x) 

(29) 
(30) 



Degree 18 transformations between hypergeometric equations with the local exponent 
differences (1/2, 1/3, 1/7) and (1/7, 1/7, 2/7) are compositions of degree 9 and degree 
2 transformations. The intermediate hypergeometric equation has the local exponent 
differences (1/2,1/7,1/7). To get a hypergeometric identity, one can compose formula 
(3) with a = 3/14, 6=1/2 and formula (25). 

Degree 24 transformations between hypergeometric equations with the local exponent 
differences (1/2, 1/3, 1/7) and (1/7, 1/7, 1/7) are compositions of degree 8 and degree 3 
transformations. The intermediate hypergeometric equation has the local exponent dif- 
ferences (1/3, 1/3, 1/7). Note that we have here a composition of two pull-back transfor- 
mations which do not admit a Coxeter decomposition, but the composite transformation 
does admit a Coxeter decomposition. To get a hypergeometric identity, one can compose 
formula (20) with c = 2/7 and formula (23); see [Vid04, formula (76)]. 

A covering for degree 10 pull-back transformations between hypergeometric equations 
with the local exponent differences (1/2, 1/3, 1/8) and (1/3, 1/8, 1/8) is given by 



ip 4 (x) = 



Ax(x- 1) (8/3a; + 7 - 4/3) 8 



(2048/3a; 3 - 3072/3a; 2 - 3264a; 2 + 912/3a; + 3264a; + 56/3 - 17) 3 ' 



(31) 



Here /? satisfies ft 2 + 2 = 0. Hypergeometric identities are: 



2F1 



/ 5/24, 13/24 
I 7/8 



2 r l 



(5/24, 1/3 
I 2/3 



x 2 Fi 
= (l-a;)" 1 / 48 



16(4-17/3) _ 64(167-136/3) 2 , 2048(112-17/3) 3 N 



-1/16 



( 1/48, 17/48 
7/8 



('- 



8+7/3 
16 



243 

ip 4 (x) 



6561 



(32) 



1/6 



2 r l 



1/48,7/48 
2/3 



ipi{l/x) 



(33) 
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Degree 12 transformations between hypcrgcomctric equations with the local exponent 
differences (1/2, 1/3, 1/8) and (1/4, 1/8, 1/8) are compositions of a degree 3 transforma- 
tion and two quadratic transformations. The intermediate hypergeometric equations have 
the local exponent differences (1/2, 1/4, 1/8) and (1/2, 1/8, 1/8). To get a hypergeometric 
identity, one can compose formula (3) with a — 1/4, b = 1/2, formula (2) with a = 1/8, 
b = 5/8, and formula (18) with a = 1/16. 

Degree 12 transformations between hypergeometric equations with the local exponent 
differences (1/2, 1/3, 1/9) and (1/9, 1/9, 1/9) are compositions of degree 4 and degree 
3 transformations. The intermediate hypergeometric equation has the local exponent 
differences (1/3,1/3,1/9). To get a hypergeometric identity, one can compose formula 
(20) with c = 1/3 and formula (19) with a = 1/12. 

A covering for degree 6 pull-back transformations between hypergeometric equations 
with the local exponent differences (1/2, 1/4, 1/5) and (1/4, 1/4, 1/5) is given by 



tp 5 (x) = 
Hypergeometric identities are: 

*( 3/20 3/ rh) - o-^ 



4i x (x - 1) (4x - 2 - Hi) 4 



ix - 4 + 3i) 5 

1/40, 9/40 



jFi 



3/4 



/ 3/20, 2/5 
2 H 4/5 



= {i- x )-y^{i-^i x ) 



-1/10 



Fi 



tp 5 (x) 

1/40,11/40 
4/5 



(34) 



(35) 



(36) 
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